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We report on magneto-optical studies of Bi 2 Se 3 , a representative member of the 3D topological 
insulator family. Its electronic states in bulk are shown to be well described by a simple Dirac- 
type Hamiltonian for massive particles with only two parameters: the fundamental bandgap and 
the band velocity. In a magnetic field, this model implies a unique property - spin splitting equal 
to twice the cyclotron energy: Es — 2Ec. This explains the extensive magneto-transport studies 
concluding a fortuitous degeneracy of the spin and orbital split Landau levels in this material. The 


Eg = 2Ec match differentiates the massive Dirac 

electrodynamics, for which Eg = Ec always holds. 

PACS numbers: 71.70.Di, 76.40.+b, 78.30.-j, 73.20.-r 

Inspiring analogies to relativistic systems have largely 
helped to elucidate the electronic properties of two- 
dimensional graphene mm, surface states of topologi¬ 
cal insulators (TIs) [3H5], novel three-dimensional (3D) 
semimetals 1311 as well as certain narrow gap semicon¬ 
ductors m- Here, we report on magneto-optical studies 
of bulk Bi 2 Se 3 , which imply the approximate applica¬ 
bility of the Dirac Hamiltonian for massive relativistic 
particles to approach the band structure of this popular 
representative of the TI family. 

The dispersion relations of genuine massive Dirac 
fermions in quantum electrodynamics are defined by two 
parameters: the energy gap 2A between particles and 
antiparticles and velocity parameter vd- At low ener¬ 
gies, i.e., in the non-relativistic limit, these dispersions 
become parabolic and characterized by the same effec¬ 
tive mass rriD = (rest Dirac mass). Such disper¬ 

sions resemble the cartoon sketch of a direct gap semi¬ 
conductor, which may be conventionally described using 
Schrodinger equation, completed by extra Pauli terms in 
order to include the spin degree of freedom. In contrast, 
no additional terms are needed when Dirac equation is 
employed, since it inherently accounts for spin-related 
effects. For instance, when the magnetic field B is ap¬ 
plied, Dirac equation describes both cyclotron {Ec) as 
well as spin {Eg) splitting of the electronic states and 
implies that these two splitting energies are the same 
and linear with B in the non-relativistic approximation: 
Eg = Ec = hujc = heBlrriD = ehBv^/A. For free elec¬ 
trons, this condition is equivalent to the effective g factor 
of 2 in Eg = gjiiBB (Bohr magneton gs = eh/2nio) pH] . 

In this Letter, we demonstrate experimentally that the 
conduction and valence bands of Bi 2 Se 3 are both, with 


electrons in bulk Bi 2 Se 3 from those in quantum 


a good precision, parabolic (perpendicular to the c-axis) 
and characterized by approximately the same effective 
mass. This crucial observation implies a great simplifi¬ 
cation of the multi-parameter Dirac Hamiltonian um 
commonly used to describe the bands of this material. 

The resulting simplified Dirac Hamiltonian differs from 
that of the genuine quantum electrodynamics system 
only by relevant (additional) diagonal dispersive terms, 
and importantly, it remains to be defined by two parame¬ 
ters only: by the bandgap energy 2A and velocity param¬ 
eter vd- These are directly read from our optical exper¬ 
iments, or alternatively, the vd parameter may be taken 
from the measurements of the Bi 2 Se 3 Dirac-cone surface 
states m- Remarkably and in contrast to genuine Dirac 
fermions, the electrons in Bi 2 Se 3 approximately follow 
the rule that their spin splitting is twice the cyclotron 
energy Eg = heBvj^/A = 2huJc = 2Ec. The effective 
mass, common for carriers in the conduction and valence 
bands, is thus roughly rric = my = 2A/u|) = 2m d and 
the spin splitting expressed in terms of the effective g 
factor, ge = Oh = 2m{^jmB‘ Our simplified view of the 
bands in Bi 2 Se 3 is not perfect (departures are extensively 
discussed), though it accounts well for the present ex¬ 
perimental results as well as for a number of magneto¬ 
transport data reported in the past and agrees with the 
recent estimate of the electron g factor. 

The presented experiments have been performed on a 
290-nm-thick layer of Bi 2 Se 3 grown by molecular beam 
epitaxy on a semi-insulating InP(lll)B substrate 
for data obtained on another specimen prepared under 
analogous conditions see Supplementary materials na. 
The 3D structure of the InP surface (with 2 nm root 
mean square roughness) transfers the stacking order of 
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Fig. 1. (color online) Far infrared magneto-transmission 
spectra of the Bi 2 Se 3 specimen. The infrared active phonon 
modes a and /3 are at higher magnetic helds accompanied by 
CR absorption, which follows linear in B dependence, see the 
inset, and implies the effective mass of electrons rUe = (0.140± 
0.005)mo. The plotted transmission Tb was normalized by 
that of the bare InP substrate. 

the substrate to the epilayer, resulting in a complete sup¬ 
pression of twinning, which is otherwise a ubiquitous de¬ 
fect in Bi 2 Se 3 thin films. The after-growth annealing in 
the Se atmosphere, reducing the final density of Se vacan¬ 
cies, helped to keep the electron density below 10^^ cm“^ 
(with the mobility /i in the 10^ cm^.V“^.s“^range), as 
confirmed in magneto-transport experiments. Impor¬ 
tantly, the thin epitaxial layer of Bi 2 Se 3 enabled trans¬ 
mission experiments at photon energies above the funda¬ 
mental interband absorption edge of this material. 

To measure the magneto-transmission spectra, a macro¬ 
scopic area of the sample (^4 mm^) was exposed to the 
radiation of a globar, which was analysed by a Fourier 
transform spectrometer and, using light-pipe optics, de¬ 
livered to the sample placed in a superconducting or re¬ 
sistive magnet. The transmitted light was detected by 
a composite bolometer placed directly below the sample, 
kept at a temperature of 1.6 K. All measurements were 
done in the Faraday configuration with light propagating 
along the c axis of Bi 2 Se 3 {z axis). In experiments per¬ 
formed with circularly polarized light, a glass linear po¬ 
larizer and a zero-order MgF 2 quarter wave plates (cen¬ 
tered at A = 4 or 5 /im) were used. 

The optical response of the thin Bi 2 Se 3 layer has been 
probed in both far and middle infrared spectral regions. 
At low energies, the response, see Fig. is dominated 


by infrared active phonon modes a and /3, which ex¬ 
hibit a weak coupling to the magnetic field [T7]. At 
higher magnetic fields, cyclotron resonance (CR) absorp¬ 
tion is well formed and it disperses linearly with B. The 
slope of this dependence provides us with an estimate of 
the electron effective mass: rrie = (0.140 ± 0.005)mo, 
which well falls into a relatively broad range of val¬ 
ues, rrie = (0.12 — 0.16)mo, deduced from other ex¬ 
periments [TSH22]- The interband absorption of Bi 2 Se 3 
exhibits a fairly rich response in magnetic fields, see 
Fig-H a). Firstly, at low 5, a distortion of the absorp¬ 
tion edge at the energy slightly above 200 meV appears. 
At higher fields {B > 10 T), the quantum regime is 
approached {fi.B > 1) and a series of interband inter- 
Landau level (inter-LL) resonances emerges. These reso¬ 
nances are almost equidistant in energy and follow nearly 
linear in B dependence, see Fig. [^b). 

The observed linearity of the optical response in 5, in 
reference to intraband (CR absorption) as well as inter¬ 
band inter-LL excitations, points towards parabolic pro¬ 
files of both conduction and valence bands. Let us rec¬ 
oncile this crucial experimental fact with the standard 
theoretical model of electronic bands in TIs from Bi 2 Se 3 
family Ena. Using a basis of spin-degenerate Se- and 
Bi-like p-orbitals, the authors of Refs. 13 na propose a 
3D Dirac Hamiltonian (4x4), expanded to include the 
electron-hole asymmetry, uniaxial anisotropy (along the 
c axis), and importantly, the band inversion, giving thus 
rise to the TI phase (via dispersive diagonal terms). 
Since these are the kz = ^ states, which provide the 
dominant contribution to the magneto-optical response 
studied in our experiments, the situation further sim¬ 
plifies. The 3D Dirac Hamiltonian decouples into two 
complex-conjugate 2D Dirac-type Hamiltonians h and /z* 
written in the basis of |Se |), |Bi t) and |Se t), |Bi |), re¬ 
spectively: 

A + (C + M)fc2 nvDk+ \ 

hvDk- -A + (C - M)fc2 ) ’ 

where k± = kx^ iky. This Hamiltonian implies, in gen¬ 
eral, non-parabolic conduction and valence band profiles: 

^(A:) = ± ^(A + Mfc2)2 + rflv%k^, (2) 

each exhibiting up to three local extremal points, de¬ 
pending on the strength of the interband coupling vd 
(effective speed of light), the electron-hole asymmetry 
parameter C, and the diagonal dispersive term M (neg¬ 
ative for systems with the band inversion). The basis 
of the h Hamiltonian allows us to associate a given spin 
projection to each band: Ed and £y^. To satisfy the 
time-reversal and inversion symmetries of Bi 2 Se 3 , the 
Hamiltonian h* provides an analogous solution, with the 
spin projections rotated, Ed and so we finally obtain 
twice spin-degenerate conduction and valence bands. 
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Fig. 2. (color online) Part (a): Rela¬ 
tive transmission spectra of Bi 2 Se 3 in the 
middle infrared spectral range plotted for 
selected values of R. At R = 32 T, in¬ 
dividual excitations are denoted by ver¬ 
tical arrows and identified by the corre¬ 
sponding index n. Part (b): Positions of 
experimentally observed interband exci¬ 
tations as a function of B. The dashed 
lines represent theoretical ht to data de¬ 
scribed in the text. The inset shows the 
magneto-transmission spectrum taken at 
R = 25 T measured with a defined cir¬ 
cular polarization of light. Notably, it is 
the normalization by To, which induces 
the modulation of Tg/To curves around 
the zero-field interband absorption edge 
oihuj = 2A-\-EF{l-\-me/mh) ~ 225 meV. 


To make the above dispersions (§ parabolic in a broad 
range of energies, z.e., to make the model consistent 
with our magneto-optical data, the specific condition 
= —4MA has to be satisfied. Notably, this is 
only possible for systems in the TI phase when M < 0 
(by definition A > 0). The bands then take a simple 
form, Tc = A + (C — M)k‘^ and Sy = —A + (C + M)k‘^^ 
and are characterized by well-defined effective masses: 
rrie = /[2{C — M)] = 21i?/{h?/rriD + 4(7) and rrih = 

—h‘^/[2{C + M)] = 2/i^/ {li?/mjj — 4(7), for electrons and 
holes, respectively. Interestingly, the corresponding re¬ 
duced mass equals to the Dirac mass: l/rue + = 

XjvfiD = v]jIA. Clearly, in the case of a relatively weak 
electron-hole asymmetry [C |^|), the expressions fur¬ 
ther reduce to rUe ~ rrih ~ 2m d- 

When the magnetic field is applied, the bands in Bi 2 Se 3 
transform into Landau levels (LLs). The Dirac-type 
Hamiltonians h and /i* give rise to particular electron 
and hole zero-mode LLs: To,e = A + ((7 + M)eB/h and 
Eq ^ = —A + ((7 — M)eB/h. These zero-mode levels 
are typical of TIs (see, e.g., Refs. 0) - they disperse 
strictly linearly with H, they are spin polarized, insensi¬ 
tive to the strength of the interband coupling vd and 
they cross each other at the field of By = hA/\eM\. 
The LLs with higher indices {n > 0) follow, assum¬ 
ing parabolic bands with a relatively weak electron-hole 
asymmetry (7 <C \M\ (a posteriori justified by our ex¬ 
perimental data), nearly linear in B dependence. For h 
and h* Hamiltonians we get the LL spectrum: = 

Eo,e + T* g = Tn-i,e + , E* j^ = Eq^^ - 

and En^h = ^n-i.h ~ ^ where = eBlme,h,D- 

Importantly, the shift huj^ between the two LL series 
corresponds to the spin splitting Eg = heBvj^/A, which 
may be expressed also in terms of a ^ factor: Qe = = 

2mo/mg), see Supplementary materials [15]. In analogy 


to massive particles in quantum electrodynamics, this 
spin splitting is given just by the energy bandgap (2A) 
and the effective velocity of light (vd) and it is the same 
for electrons and holes (particles and antiparticles). On 
the other hand, the effective masses of electrons and holes 
depend on the diagonal terms M and (7 and this implies 
a certain ratio, EsjEy > 1, between the cyclotron en¬ 
ergy and spin splitting (notably, Ec/Eg = 1 always holds 
for free electrons in vacuum). The spin-splitting of elec¬ 
tronic bands in TIs (with M < 0) should thus manifest 
at lower magnetic fields, prior to Landau level quantiza¬ 
tion cf. Figs.|^a-c). This, in fact, accounts for the initial 
distortion of the absorption edge observed at low mag¬ 
netic fields, see Fig. [^a). Interestingly, for rather small 
electron-hole asymmetry ((7 ^ \M\ ^ me 2m d)^ we 
get 2Ee = 2huJc = heBjmD = geEsB = Eg. 

The Dirac Hamiltonian 0 in a magnetic field gives 
rise to two series of dipole-active inter-LL excitations 
n^n 1 and n^n — 1, active in and (t~ polar¬ 
ized light, respectively. This has been corroborated both 
experimentally and theoretically, for instance, in the con¬ 
text of graphene [23l[23]. However, a closer look at the 
corresponding matrix elements shows that, in a gapped 
system (A ^ 0), n^n — 1 transitions dominate inter¬ 
band inter-LL absorption within the spectrum of the h 
Hamiltonian, and vice versa, the n ^ n + 1 series in the 
h* Hamiltonian, see Fig. |^c) for illustration and Sup¬ 
plementary materials [T5| for details. This behavior may 
be viewed as a spin-dependent optical activity, implying 
necessity to reverse spin during the interband absorp¬ 
tion. The interband excitations, Sy^ Bcf 
connecting states within /i* or h Hamiltonians, respec¬ 
tively, are thus active in and (j~ polarized radia¬ 
tion only. Notably, a similar situation is encountered in 
gapped graphene and transition-metal dichalcogenides. 
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Fig. 3. (color online) Part (a): Approximate profile of electronic bands in Bi 2 Se 3 at the center of the Brillouin zone (for kz = 0). 
The dashed lines show the (dispersive) diagonal term of the Hamiltonian Q. Part (b): Spin splitting of the conduction and 
valence bands {ge ~ Qh) at B = 8 T (LLs not yet resolved g.B < 1). Part ^): Fan chart of nearly linear in B LLs in Bi 2 Se 3 . 
The approximate match of LLs En+i,e ~ L^n,e ^ En^h ~ is a direct consequence of the Dirac-type model for electronic 

states in Bi 2 Se 3 discussed in the text. Vertical arrows show selected electric-dipole-active interband inter-LL resonances, CR 
absorption as well as (expected) magnetic-dipole-active electron spin resonance (ESR) absorption in the quantum limit. Part 
(d): High-magnetic-field extrapolation of zero-mode LLs, with the crossing point at Be ~ 300 T. 


where valley-sensitive selection rules for circularly polar¬ 
ized light appear [25]. Let us also note that intraband 
(CR) absorption is always active in and a~ polarized 
radiation for electrons and holes, respectively. 

The LL spectrum and the selection rules allows us to 
identify individual resonances in the interband response 
in Fig.|^ When the quantum limit is reached {B > 20 T), 
the lowest in energy observed absorption line is the 
Lq -^Ll transition, active in cr+ polarized light. The 
parent Li ^Lq line, active in a~ polarized radiation, 
does not appear, since the bottom of Eq level is always 
occupied in the n doped system, cf. the inset of Fig.j^b) 
and Fig.j^c). At higher energies, we get a series of tran¬ 
sitions, L* {n > 0), which are 

for a given n nearly degenerate in energy and active in 
cr+ and (T~ polarized light, respectively. Their spacing, 
h{uj^ + ^ allows us to read the Dirac mass 

directly from the data: ttid = (0.080 ± 0.005)mo. No¬ 
tably, this Dirac mass implies ge = = 2mQ/mD ~ 25, 

which is in very good agreement with the value gf^^ = 
27.5 [26], derived recently using spin resonance measure¬ 
ments. The Dirac mass mjj, together with the bandgap 
2A = (200 ± 5) meV, read from the low-field extrapola¬ 
tion of interband inter-LL resonances in Fig.j^b), imply 
the velocity vd = = 4.8 x 10^ m/s, in perfect 

agreement with majority of ARPES studies [13127]. 

Comparing the estimated Dirac mass with the electron 
mass (as deduced from CR absorption), we conclude that 


nie ~ 2mD- This indicates rather weak electron-hole 
asymmetry in Bi 2 Se 3 {C C \M\ = h^vj^/AA). Ne¬ 
glecting this asymmetry completely, we remain with the 
Hamiltonian 0 with two independent parameters only: 
A and vd- Interestingly, these two parameters, which can 
be easily read from infrared transmission experiments, 
fully describe the band structure of Bi 2 Se 3 : the energy 
bandgap of 2A, effective masses rrie ~ rrih ~ 2m d = 
2A/'z;|), as well as g factors, ge = gh = 2mo/mD = 
2movjylA. Notably, this match between twice cyclotron 
energy and spin splitting {2heBlme = gel^sB) has been 
in the past found as a purely empirical fact in quantum 
oscillation experiments on Bi 2 Se 3 , see, e.g.^ Refs. [^ 1^. 
Here we show that this surprising match is not acciden¬ 
tal and it straightforwardly follows from the Dirac-type 
Hamiltonian Q applied to TI with a weak electron-hole 
asymmetry and nearly parabolic bands. 

It should be also mentioned that within our “parabolic 
view” of electronic bands, Bi 2 Se 3 clearly becomes a 
direct-gap semiconductor. This is in agreement with re¬ 
cent experimental (ARPES) and theoretical studies, see, 
e.g.^ Refs. |3Qll33] « however, in contradiction with other 
ARPES data, see, e.^.. Refs. ns 127], in which the ob¬ 
served camel-back profile of the valence band indicated 
an indirect band gap. Intriguingly, our experiments, to¬ 
gether with other optical studies performed on bulk or 
thin-film specimens, see, e.g.^ Refs. |34H36] . provide a 
significantly lower band gap (200 meV) as compared to 
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values deduced from ARPES experiments (^300 meV), 
see Supplementary materials m for further discussions. 
Now we will discuss the limits of our simplified two- 
parameter model (with parabolic bands and full electron- 
hole symmetry), by confronting it with more detailed 
analysis of our experimental data. The real band struc¬ 
ture of Bi 2 Se 3 may deviate by (i) appearance of the 
electron-hole asymmetry and (ii) the departure of bands 
from exact parabolicity. The electron-hole asymmetry is 
clearly demonstrated by rrie < rrih {me < 2m which 
translates into C = (3 ± 0.5) eV.A^, but also by the dif¬ 
ference in the corresponding g factors. The latter may 
be read from a small, but noticeable, splitting between 
L* Ln+i transitions in the spectra 

taken with a defined circular polarization of light, see 
the inset of Fig. [^b). It implies ge — gn ^ 3, which may 
be explained as the contribution of the free-electron Zee- 
man term and influence of remote bands, described by, 
e.^., the Roth’s formula [HIST]. 

The deviations from bands’ parabolicity imply the depar¬ 
ture of inter-LL resonances from their linearity in B. In¬ 
deed, the transitions at higher energies and/or for higher 
LL indices, see Fig. |^b), slow down to a weak sublin- 
ear dependence. To describe this behavior, we have used 
the full (non-linearized) expressions for LLs, see Supple¬ 
mentary materials [T3], to fit the positions of individual 
resonances. We varied parameters M and A, while 
fixing C = 0, which has rather weak impact on the in¬ 
terband response. The best agreement is obtained for 
VD = (0.47 ± 0.02) X 10^ m.s-\ A = (0.100 ± 0.002) eV 
and M = —(22.5 ± 1.0) eV.A^. We may thus conclude 
that the condition = —4MA is fulfilled within a 

few percent, which validates our view of parabolic bands 
in Bi 2 Se 3 . Moreover, since C/\M\ ^ 1/10, the system 
indeed exhibits rather high electron-hole symmetry. 

The deduced strength of dispersive diagonal terms, M 
and A, allows us to estimate the critical field Bc^ at which 
the zero-mode LLs cross each other, see Fig. [^d). At 
this magnetic field, Bi 2 Se 3 changes into a semi-metallic 
(= gapless) material, for which the extended 3D Dirac 
Hamiltonian implies, see Supplementary materials m, 
the linear in kz bands, E{kz) = -EhvD\kz\^ with a high 
(LL) degeneracy = eB jh. The velocity vd is supposed 
to be slightly lower as compared to vd due to the uniax¬ 
ial anisotropy of Bi 2 Se 3 HHiiS]. The pretty high value 
of Be ~ 300 T, at the limit of currently available (semi¬ 
destructive) pulsed-field techniques [38] , makes the explo¬ 
ration of this interesting critical point difficult. However, 
this crossing field Be is expected to be lower in other TIs 
from the Bi 2 Se 3 family, with a lower band gap, e.^., in 
Bii_a,Sba, for rather low Sb concentrations [3|. 

In conclusion, we have shown that the band structure 
of Bi 2 Se 3 can be, in very good approximation, described 
by a simple Dirac-type Hamiltonian with only two free 
parameters: the effective velocity parameter vd and the 
band gap 2A. This simplified model provides us with 


reasonable estimates for both effective masses (rUg ~ 
mh ~ 2A/u|)) and corresponding g factors {ge ^ gh ^ 
2mou|)/A), and implies, for charge carriers in Bi 2 Se 3 , a 
surprising match between the cyclotron energy and spin- 
splitting: Es ~ 2Ec. Notably, this relation has been de¬ 
duced from quantum oscillations experiments performed 
on Bi 2 Se 3 in the past, but only as a purely empirical fact. 
Here we show that this directly follows from the Dirac- 
type Hamiltonian applied to a TI with nearly parabolic 
bands and a high electron-hole symmetry. 

The work has been supported by the ERC-AG projects 
MOMB and 3-TOP. Authors acknowledge discussions 
with D. M. Basko, O. Ly and M. O. Goerbig. T. B. ac¬ 
knowledges the support from the Austrian Science Fund 
(FWF), Grant No. M 1603-N27. E. M. H. and and C. 
M. thank DFG grant HA 5893/4-1 within SPP 1666. 


* milan.orlita@lncmi.cnrs.fr 

[1] K. S. Novoselov et al., Nature 438, 197 (2005). 

[2] Y. B. Zhang et a/., Nature 438, 201 (2005). 

[3] B. A. Bernevig and S.-C. Zhang, Phys. Rev. Lett. 96, 
106802 (2006); M. Konig et al, Science 318, 766 (2007). 

[4] D. Hsieh et al, Nature 452, 970 (2008). 

[5] H. Zhang et a/.. Nature Phys. 5, 438 (2009). 

[6] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010): X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 
1057 (2011) 

[7] Z. K. Liu et al., Science 343, 864 (2014). 

[8] M. Orlita et al., Nature Phys. 10, 233 (2014). 

[9] Z. K. Liu et al.. Nature Mater. 13, 677 (2014). 

[10] W. Zawadzki, Adv. Phys. 23, 435 (1974). 

[11] V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, 
“A course of theoretical physics, vol. 3, p. 455,” (Perga- 
mon, 1977). 

[12] C.-X. Liu et al, Phys. Rev. B 82, 045122 (2010) 

[13] Y. Xia et al.. Nature Phys. 5, 398 (2009). 

[14] N. V. Tarakina et al., Adv. Mater. Int. 1, 1400134 (2014) 

[15] See Supplemental Materials. 

[16] O. B. O. Ly, Eleetron spin resonanee in topologieal insu¬ 
lators: Theoretieal study. Master’s thesis, Universite de 
Strasbourg (2014). 

[17] A. D. LaForge et al, Phys. Rev. B 81, 125120 (2010) 

[18] H. Kohler and E. Wochner, phys. stat. sol. (b) 67, 665 
(1975); K. Eto et al, Phys. Rev. B 81, 195309 (2010) 

[19] J. G. Analytis et al, Phys. Rev. B 81, 205407 (2010) 

[20] N. P. Butch et al, Phys. Rev. B 81, 241301 (2010). 

[21] A. N. Sushkov et al, Phys. Rev. B 82, 125110 (2010). 

[22] H. Cao et al, Phys. Rev. Lett. 108, 216803 (2012) 

[23] M. L. Sadowski et al, Phys. Rev. Lett. 97, 266405 (2006) 

[24] V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. 
Rev. Lett. 98, 157402 (2007). 

[25] W. Yao, D. Xiao, and Q. Niu, Phys. Rev. B 77, 235406 
(2008); F. Rose, M. Goerbig, and F. Piechon, ibid. 88, 
125438 (2013) 

[26] A. Wolos et al, AIP Conf. Proc. 1566, 197 (2013). 

[27] Z.-H. Zhu et al, Phys. Rev. Lett. 107, 186405 (2011) 

[28] H. Kohler and H. Fischer, phys. stat. sol. (b) 69, 349 
(1975) 



6 


[29] B. Fauque et al, Phys. Rev. B 87, 035133 (2013) 

[30] Y. L. Chen et al, Science 329, 659 (2010) 

[31] O. V. Yazyev et al, Phys. Rev. B 85, 161101 (2012) 

[32] I. Aguilera et al, Phys. Rev. B 88, 045206 (2013) 

[33] I. A. Nechaev et al, Phys. Rev. B 87, 121111 (2013) 

[34] D. Greenaway and G. Harbeke, J. Phys. Chem. Solids 
26, 1585 (1965). 


[35] H. Kohler and J. Hartmann, phys. stat. sol. (b) 63, 171 
(1974). 

[36] K. W. Post et al, Phys. Rev. B 88, 075121 (2013), 

[37] L. M. Roth, B. Lax, and S. Zwerdling, Phys. Rev. 114, 
90 (1959) 

[38] O. Portugall et aL, C. R. Phys. 14, 115 (2013), and 

references therein. 


Supplementary Information for 


Magneto-optics of massive Dirac fermions in bulk Bi2Se3 

by M. Orlita, B. A. Piot, G. Martinez, N. K. Sampath Kumar, C. Faugeras, M. Potemski, 
C. Michel, E. M. Hankiewicz, T. Brauner, C. Drasar, 

S. Schreyeck, S. Grauer, K. Brunner, G. Gould, G. Briine, and L. W. Molenkamp 


In this supplementary material, we present details of Landau level (LL) spectrum and optical selection rules for bulk 

Bi2Se3. 


BULK LANDAU LEVELS 


Liu et al. proposed a 3D Dirac Hamiltonian to describe the bulk states in Bi 2 Se 3 [S 1 ,S 2 ]. This Hamiltonian, written 
in the basis {|Se A 5 |Bi t) JSe t), |Bi A}? mads: 


H — e(/s)l4x4 + 


/ M{k) 
A{k)k- 
0 

\-B{k,)k, 


A{k)k-^ 0 

-M{k) B{k,)k, 
B{k,)k, M{k) 
0 A{k)k-^ 


-B{k,)kA 

0 

A{k)k- 

-M{k) J 


(3) 


where €{k) = Cq + Cik^ + M{k) = Mq + Mi+ M 2 k^, B{k,) = Pq + B 2 kl, A{k) = Aq + A 2 k^, k± = k^± iky 

and k‘^ = k‘^ -\- ky. In the simplified model, employed in the main text, we restrict this Hamiltonian by neglecting the 
k^ terms and redefine A = Mq, M = M 2 , C = C 2 , Go = 0 and vd = Aq/H. 

To examine the Landau level structure, we introduce the magnetic field (along the 2 : axis) by means of the Peierls 
substitution k —tt := /c + where A is the vector potential, B = V x A and B\\z. In the quantizing magnetic 
fields, the LL spectrum of electrons and holes takes the form of En^e{kz) and En^h{kz)- The inversion symmetry of 
the system implies: E^^^^jy^ikA) = Bn^elh{~kz) and therefore dE^^elhl^kz = 0 ai kz = 0. This induces a series of 
singularities in the density of states for electrons and holes, but importantly for our case, also in the joint density 
of states (for each LL and inter-LL resonance, respectively). Using a parabolic expansion in kz (approximately valid 
for kz ^ 0)^ we get the characteristic pe-/i(^) 1/A?iuj — E profiles in the vicinity of each resonance in the joint 
density of states. Such a profile is typical of ID systems with a parabolic dispersion. The singularities in pe-h^ m 
reality smoothed by disorder, give rise to the experimentally observed inter-LL resonances. In principle, more complex 
Landau level profiles in the momentum kz may provide further singularities/resonances (due to states at kz 7 ^ 0), 
nevertheless, we did not identify any such transitions in our magneto-transmission spectra. The magneto-optical 
response of Bi 2 Se 3 {B\\z) is thus dominantly determined by = 0 states. 

For kz = the 3D Hamiltonian © reduces to two 2D Dirac-like Hamiltonians: 


Ho = 


(hoin) 

V 0 


0 ) 
h^{n)) 


with ho{iT) = e( 7 r)l 2 x 2 


fM{Tz) AqTT+ \ 

\Ao7r_ 


(4) 


that will be referred to as the “full model” in the following discussion. Now, without kz terms, we have e( 7 r) = 
where C breaks the particle-hole symmetry. M{n) is the k dependent mass term M{n) = A -f- with A 

determining the band gap Eg = 2 A at k = 0. 

To study the formation of LLs, we introduce the ladder operators: 




( 5 ) 
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with the magnetic length Ib = y These operators obey the standard relations [a, a^] = 1, a \n) = \n — 1) and 
a’t’ \n) = y/n + 1 |n + 1). Using the definitions (§, we can rewrite ho{Tr) in terms of the raising and lowering operators: 




-{C + M){a^a- 




ip-oa' 

1{C -M){aU+\) 


(6) 


The form of this Hamiltonian suggests the following ansatz for the eigenstates: 

= (c„2ln-\)) ^ (^0^) ^ 

Solving the Schrodinger equation, we find the LL spectrum: 


M 

En ,. c\ — T 2 rpy Ti 

C + M 


^n,a — j2 



c , M 

W + A + 2-^n 

'■B ‘B / '■B 




Eo,e ■ 


/2 


( 8 ) 


where 5e = +1 for electrons and = —I for holes. Note that each state with the energy of ^n>i,a is always a 
superposition of the |Se, n, 1) level and the |Bi, n — 1, t) level. In contrast, the ^o,e state is fully polarized as a 
|Se, D level. The presence of such polarized zero-mode LLs is characteristic of Dirac-type Hamiltonians and it is 
well-known, e.^., from physics of graphene [S3]. However, in graphene, the zero-mode levels (at K and K' points) are 
polarized in pseudospin not in real spin as in the case of Bi 2 Se 3 . 

The Landau levels for the Hamiltonian are found in an analogous way. Here, we take 

" (| 0 )) 

as the ansatz for the eigenstates and get LLs for the Hamiltonian: 


77 '* _ 

En.a 


M 

IT +^j^n + So 


c , M V O^o 

72 

■ ^B / ^B 


771 * _ 

Eo^h — 


-A 


C-M 

P * 


( 10 ) 


Note that just denotes energies/states belonging to the Hamiltonian and does not stand for complex conjugation. 
We will keep this notation in the next sections. For the Hq Hamiltonian, the state with the energy F^n>i,a is always 
a superposition of the |Se, n — 1, t) level and the |Bi, n, 1) level. Again, the zero-mode LL is fully spin-polarized, in 
this case as a |Bi, 1) state. 

The LL spectrum of the Hamiltonian Q (calculated within the full model) is plotted in Fig. [^for parameters derived 
from our magneto-optical experiments, together with electric-dipole transitions discussed later on in detail. 

To simplify the LL spectra of the ho and Hq Hamiltonians, we assume perfectly parabolic bands {Aq = = —4MA) 

and expand Eq. (§ and Eq. ( pu] ) for small magnetic fields. This way we get LLs strictly linear in the applied magnetic 
field: 


En,e — Eq^c T huJ^Tl^ 

-E'n+l.e =En,e + 

- Swf, 

E*n,h =Elh - Hn, 

where the cyclotron frequencies are defined as ujc^^ = eB and = eB/mn with the effective masses rrig/fi = 

2h‘^/ {jP/rriB ±4(7) and the Dirac mass rriB = A/i;|>, respectively. Importantly, for parameters deduced from our 
experimental data, this simplified LL spectrum is nearly identical to that calculated within the full model, see Eig. 
Eurther information about Landau levels may be obtained from the analysis of individual eigenstates, which (with the 
exception of zero-mode levels) represent a superposition of spin-up and spin-down states in the conduction and valence 









Fig. 4. The Landau level spectrum of the full model Blue levels originate from the ho sub-Hamiltonian, red levels from 
the ho Hamiltonian. The arrows show the dominant type of electric-dipole transitions for the respective Hamiltonian. 




B (T) B (T) 


Fig. 5. LLs calculated within the full model, Eq. 
for small magnetic field (z.e., by linearization in B 


see Eq. 

and high LL indices. Eor simplicity, only the levels of the ho Hamiltonian are shown. 


and (|10|), compared to LLs obtained by expansion of the same formulas 
0- Significant deviations appear only at high magnetic helds 


bands. Nevertheless, for low-energy/low index levels we may always find the dominant state in this superposition, 
see Fig. We may conclude that the eigenstates of the ho Hamiltonian with the energies of En^e and En^h can 
be considered as selenium-like spin-down and bismuth-like spin-up levels, respectively. Analogously, 

states with the energies of E’* ^ and E* ^ have Se-like spin-up {E^^\) and bismuth-like spin-down character, 

respectively. This assignment of spin-projections will facilitate the definition of g factors in the next section. 


DEFINITION OF G FACTORS 

For B = OT the energy states in our system are spin degenerate as required by time- and inversion- symmetry of 
Bi 2 Se 3 . When the magnetic field is applied, this spin degeneracy of states is lifted due to the magnetic moment fjL 
of electrons. This splitting may be described by a corresponding effective g factor in the Zeeman term, Ez = 
with fi = g^Bs/h. The total g factor comprises three contributions. As shown in our experiments, and by the 
subsequent data analysis, the main contribution results from the strong the spin-orbit coupling in Bi 2 Se 3 , which is 
inherently included within the Hamiltonian (§. Further (minor) corrections come from the free electron go ^ 2 factor 
(free-electron Zeeman term) and a perturbative contribution from remote energy bands [S4]. 

From the previous section, we know that, for small level indices n and low energies, LLs are nearly spin polarized. 
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Fig. 6. Probability density of the energy states being Bi-like (dark color) or Se-like (light color) with corresponding spin- 
polarization. ne and nh denote the LL with index n for electron (e) or hole (h) states. This probability is shown in the part 
(a) for LLs belonging to the ho Hamiltonian, where these levels are a mixture of |Se, 1) and |Bi, t) states. The part (b) shows 
this probability for the Hq Hamiltonian, where LLs are a mixture of |Se, t) and |Bi, 1) states. Apart from the zero-mode LLs, 
one can see that the spin-polarization becomes weaker with increasing magnetic held. However for magnetic helds B < 40T 
and low level indices the dominant spin polarization of the energy states stays the same as in the B —0 T limit. 


In addition, the Landau levels in the valence and conduction bands are bismuth- and selenium-like, respectively. 
Following this fact, we can express the g factors of charge carriers, in the conduction (c) and valence (v) bands, in 
terms of LLs calculated from the Hamiltonian Eq. 0: 

9e = 9e =5Se(n, B) = - E^^,)I{ixbB\ (12a) 

9v = -9h = 9 Biin,B) = {E^''^ - E^'^)/{^bB) = {En+i,h - (12b) 


This dehnition of the factors for electrons and holes is consistent with magnetic-dipole selection rules: n ^ 1, 

which interconnect states belonging to the and Hamiltonians [S5]. In addition, we may crosscheck this dehnition 
by taking genuine Dirac particles (electrons) in the vacuum. For this, we have to take Aq = he {vd = c), A = moc^ 
and M = C = 0, where c is the speed of light and mo the rest mass of a free electron. Indeed, we get = 2 as 
expected. 

The dehnition (12) implies g factors that, in general, depend on the magnetic held as well as on the LL index. However, 
within our parabolic approximation = —4MA) and for LLs linearized in 5, see Eq. ( pTj ), the spin-splitting of 

electrons and holes becomes linear in magnetic held, Es = huj^, implying thus gc = —gv = Qe = Qh = ^mo/rriD = 
2movjylA. Eor the experimentally determined Dirac mass tud = (0.080 ± 0.005)mo we get ge = g^ ^ 25. 

Assuming the parabolic approximation = —4MA), but taking the full expressions for Landau levels ^ and 

( [1Q| ), z.e., not linearized in 5, the g factors for electrons and holes slightly differ and also gain a weak magnetic-held 
and Landau-level-index dependence, see Eig. Nevertheless, this is not sufficient to account for the experimentally 
observed difference, ge — Qh ~ 3, derived from our data in the main text. This may only be explained by further 
corrections (the Zeeman term with the free-electron go = 2 factor and the inhuence of remote bands). 


ZERO-MODE LANDAU LEVELS IN HIGH MAGNETIC FIELDS 

In high magnetic helds, the spin-polarized zero-mode LLs Ee,o = A + (C + M)//^ and = —A + (C — 
approach each other and become well separated from the rest of levels. In such a case, we can describe them by an 
effective Hamiltonian: 

/ Ee,o -hvDkz\ 

\—hvDkz Eh,o ) 


-zero—mode — 


(13) 
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Fig. 7. g factor for charge carriers in the conduction and valence bands, gdn , B) a nd gv jn, B ), exhibiting a weak dependence on 
the level index n and the magnetic held strength B according to equations (12a) and (12b). (a) With n < 5 and B < 40 T the 
approximate constant value for gc shows a deviation of < 3 % from the n- and S-dependent g factor, (b) For gv this deviation 
grows to < 12%. Let us note that we use the notation gc — ge and gv = —gh- 


which can be derived from the {kz dependent) 3D Dirac Hamiltonian proposed in Refs. [S 1 ,S 2 ], in which we neglect 
the terms cubic in terms square in kz and introduce the magnetic held via Peierls substitution. Taking account 
of the experimentally observed uniaxial anisotropy of Bi 2 Se 3 [S6,S7], one can assume that vd = Bq/H < vd- The 
effective Hamiltonian (13) is equivalent to a ID Dirac-type Hamiltonian with the band gap of F^e,o — This gap 
vanishes at the crossing held Be = h/S./\eM\ 300 T, when the system changes from a (narrow gap) semiconductor 

into a gapless semimetal. The electronic bands in such a semimetal are equivalent to the ID Dirac-type channel, 
B{kz) = 4zhvD\kz\ with a strong (LL) degeneracy of states eB/h. 


SELECTION RULES AND MATRIX ELEMENTS 


To describe the response of our system to an externally applied electromagnetic held, we employ the standard (linear- 
response) Kubo-Greenwood formalism. For polarized radiation, the optical conductivity tensor, in a system with 
eigenstates |T^) and corresponding energies reads: 


cr±(co’, B) (X i 


B 

UJ 


E 

n,n' 


^(/n - fn') 


\{^n'\v±{B)\^n)f \ 

En - En' - fu^ + ij j 


(14) 


where fn is the occupation factor, 7 the phenomenological broadening parameter and v± the velocity operators. The 
matrix elements |h±(H)|Tn) determine the active electric-dipole transitions (selection rules) between different 
eigenstates. 

The velocity operator can be directly obtained from Eq. (|^ by calculating Vi = ^§ 7 ^, where i = x,y and v± = 
{vx±iVy)jV2: 


0 0 0 

dvD 0 0 

0 0 V2 

Vo 0 "0 


( 15 ) 


Notably, these velocity operators (with two independent 2 x 2 diagonal blocks) imply that the electric-dipole transitions 
are not active between pairs of LLs belonging to different Hamiltonians ho and h^. This is different from magnetic- 
dipole transitions, which connect states originating in different Hamiltonians ho and h^. 

Taking the eigenstates |T^), z.e., LLs expressed by Eqs. 0 and ^ arranged as: 


l^n) = 



(16) 
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we may calculate the selection rules sensitive to the circular polarization of the infrared radiation: 

n' I ^+(-^) l^n) ~ ^n^n,n' — 1 ^ri^n,n' — 1 with U ^ 0, fl 7) 

{^n'\V-{B)\^ri) = + + l ^ith n > 1, 

where A, M, C, B) and T* = T*{Aq^ A, M, C, B) are amplitudes belonging to transitions within the ho 

and ho Hamiltonians, respectively. Importantly, we get the same selection rules for electric-dipole transitions between 
LLs belonging to the ho and h^ Hamiltonians, n nil. Nevertheless, the corresponding amplitudes and i** 
may strongly differ. Here we should again recall that the symbol refers to the given sub-Hamiltonian h^ and 
does not denote the complex conjugation {\J-^\ 7 ^ iJ'nl)- The difference in amplitudes is clearly seen in the matrix 
elements for interband inter-LL absorption, see Fig. The interband inter-LL absorption between levels belonging 
to the Hamiltonians ho and h^ is dominantly active in a~ and cr+ polarized radiation, respectively. 



Fig. 8. Matrix elements Fnm,± = h^ KTm| v±{B) |^n)|^ for interband inter-LL transitions in Bi 2 Se 3 plotted in the logarithmic 
scale (for m, n = 1, 2). For a~ and polarized light, interband absorption is dominated by transition between LLs originating 
from the ho and ho sub-Hamiltonian, respectively. 


ENERGY BAND GAP OF BI 2 SE 3 

The magneto-transmission experiment, presented in this paper, provides a fairly precise estimate of the energy band 
gap in Bi 2 Se 3 : 2A = (200 ± 4) meV. This result is in very good agreement with other optical studies. For instance, 
it matches well the value of ^ 175 and 160 meV expected for the band gap at low temperatures, as extracted from 
extensive reflectivity measurements on a series of bulk specimens with different electron densities in Refs. [S 8 ] and [S9], 
respectively. Similarly, our results correspond very well to conclusions of recent low-temperature infrared transmission 
studies performed on thin layers of Bi 2 Se 3 (a series of samples with thicknesses below 100 nm) prepared by molecular 
beam epitaxy on a (111) oriented silicon substrate [SIO]. 

On the other hand, ARPES studies of bulk Bi 2 Se 3 report gap values close to 300 meV, see, e.^., Refs.[Sll-13], which 
are significantly higher as compared to our results and other optical studies, which may invoke questions about the 
nature of the thin Bi 2 Se 3 samples studied in these works. It is, for instance, not a priori clear whether and how the 
substrate properties and particular growth conditions influence the observed energy band gap and the overall band 
structure. 

To provide another independent verification of the band gap value deduced optically, we have performed low- 
temperature infrared transmission measurements on thin self-standing Bi 2 Se 3 layers, prepared simply by slicing bulk 
crystals. A typical infrared transmission spectrum is plotted in Fig. This spectrum has been measured on a 
10 /rm-thick specimen, with the electron concentration close to 10^^ cm“^. This density has been deduced from 
magneto-transport measurements and it is comparable (slightly higher) with respect to the electron density in the 
thin film of 812803 studied in our magneto-transmission experiments. 

The presented transmission spectrum exhibits a fairly sharp high energy cut-off at Eg ^ 250 meV. This provides us 
with a well-defined upper hound for the electronic band gap of 812803 , 2 A < E*, as schematically depicted in the 
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Fig. 9. Low temperature infrared transmission spectrum of a 10-micron-thick free-standing layer of Bi 2 Se 3 . The inset schemat¬ 
ically shows the relation between the high energy cut-off of the transmission spectrum and the electronic band gap: Eg > 2A. 
The pronounced modulation of the spectrum corresponds to the Fabry-Perot oscillations, showing rather high crystalline quality 
of the studied Bi 2 Se 3 bulk specimen. 


inset of Fig.[^ Clearly, this high-energy cut-off is significantly below the band gap of 300 meV, which is deduced from 
ARPES measurements. Instead, in our case, the band gap should approach 2A ^ Eg — 2 Ef ^ 200 — 210 meV, as 
implied by the Burstein-Moss shift in materials with high electron-hole symmetry (rUe ^ TUh)- The Fermi level has 
been estimated as = 20 — 25 meV for the given electron density. 

To conclude, the ARPES data indicate the band gap, which is significantly higher as compared to rather direct 
optical measurements presented in this paper as well as those performed by other groups. At present, we do not have 
a clear explanation for this intriguing difference, nevertheless, we speculate that ARPES is, as a matter of fact, a 
surface-sensitive technique. As such, the deduced band gap might be influenced by specific band-bending effects on 
the samples’ surfaces, notably in the system with an inversed order of electronic bands. This difference clearly shows 
that the consensus about the size of the band gap in Bi 2 Se 3 has not yet been established. This includes also on-going 
discussions, one versus another ARPES data, about the direct/indirect nature of the band gap in this material, see 
Refs. [Sll-15]. 


ADDITIONAL EXPERIMENTAL DATA 


Here we present complementary experimental data, obtained in high-held infrared magneto-transmission experiments 
performed on a 102 nm-thick Bi 2 Se 3 layer on a InP(lll)B substrate. This sample was prepared using MBE technique 
under conditions analogous to the 290-nm-thick sample described in the main part of the paper and it is weakly 
n-doped with the electron density slightly below 10^^ cm“^. In spite of a lower signal-to-noise ratio obtained on this 
thinner sample, the observed magneto-optical response, see Eigs. pR]and El allows US to draw the same conclusions 
about the electronic band structure of Bi 2 Se 3 as in the case of the 290-nm-thick specimen. 

The overall linear in B optical response, including intraband and interband inter-LL resonances, points towards 
parabolic prohles of both conduction and valence bands. The Dirac mass and the band gap (derived from the 
separation and low-magnetic-held extrapolation of resonances in Eig. respectively) as well as the effective mass of 
electrons (read from CR absorption in Eig. 11) are nearly identical to values obtained from the 290-nm-thick sample: 
'rriD = (0.080 ± 0.005)mo, 2A = 190 ± 5 meV and rUe = (0.14 ± 0.1)mo. The parabolic prohles of electronic bands 
together with the condition rrie ~ 2m d thus imply also for this Bi 2 Se 3 sample the specihc match between spin-splitting 
and cyclotron energy {Eg = 2Ec). 

In our deeper analysis, we compared (htted) the experimentally read positions of resonances with theoretically expected 
transition energies calculated using full (non-linearized) expressions for LLs Eqs. and ( [Tq| ), see the inset of Eig. 

The best agreement was found for parameters vd = (0.45 ± 0.03) x 10^ m.s“^, A = (0.095 ± 0.003) eV and M = 
— (22.5 ± 1.5) eV.A^, which practically match those deduced from the 290-nm-thick sample, see the main text. The 
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Fig. 10. Relative transmission spectra of Bi 2 Se 3 in the middle infrared spectral range plotted for selected values of the magnetic 
field. At R = 24 and 32 T, individual inter-LL excitations are denoted by vertical arrows and identified by the corresponding 
index n. The fanchart of the observed resonances is plotted in the inset. The dashed lines represent theoretical fits described 
in the text (for C = 0). 


deviation from the condition = —4MA for the exact parabolicity of electronic bands does not exceed a few 

percent for this set of parameters. This provides us with another justification for our approximation in which we 
describe the band structure of Bi 2 Se 3 using a simplified Hamiltonian implying only two parameters: the band gap 
2A and the velocity parameter V]j. 
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Fig. 11. Far infrared magneto-transmission spectra of the 102-nm thick Bi 2 Se 3 on a InP substrate. The CR absorption is 
manifested as a dip in the relative magneto-transmission spectra Tb/Tq (denoted by vertical arrows) and follows linear in B 
dependence, which implies the effective mass of electrons rUe = (0.14 ± 0.01)mo, see the inset. At higher magnetic fields, 
the response at low energies is characterized by field-induced transmission, Tb/Tq > 1, which is due to suppressed zero-field 
Drude-type absorption, see Ref. [16] for analogous behavior in highly doped graphene. At low energies, a pronounced effect of 
interaction between CR and a and {3 phonon modes are also observed (cf. Fig. 1 in the main part of the paper). 











